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Abstract 



We study the flavour singlet mesons from first principles using lattice 
QCD. We explore the splitting between flavour singlet and non-singlet for 
^ ' vector and axial mesons as well as the more commonly studied cases of 

^f) , the scalar and pseudoscalar mesons. 

o 
o 

t^ ! 1 Introduction 

O 

The splitting in mass of flavour singlet and non-singlet mesons with the same 
quark content arises from gluonic interactions. The assumption that these are 
t^ ' small is known as the OZI rule. For the pseudoscalar mesons this splitting is 

not small (it is related to the t], tj mass difference), basically because of the 
^ i' impact of the anomaly. For scalar mesons the splitting is also expected to be 

O . large because of mixing with the nearby scalar glueball. It is usually assumed 

that the OZI rule is in good shape for the vector and axial mesons. This is 
difficult to check experimentally because of mass shifts induced by decays etc. 
Here we study this splitting from first principles using lattice QCD. We work 
7^ . in a simplified world with N f — 2 flavours of degenerate quarks. This is still 

C^ ' sufficient to explore the sign and magnitude of any splitting which can then have 

a phenomenological impact on the interpretation of the observed spectrum. In 
this N f — 2 world, isospin is exact so we classify the flavour singlet state as mo 
and the non-singlet as m\. We indeed do flnd that the splitting is largest for 
the scalar and pseudoscalar mesons, but we are able to estimate the magnitude 
of the splitting for vector and axial (both /i, a\ and /ii, h\ types of axial) 
mesons. For the vector and axial channels, the corresponding glueball masses 
are known from lattice studies |3, 0] to be much heavier and so the expectation 
is that the splitting toq — m\ due to such glueball effects would be small and 
negative (since mixed states repel). This is the same sign as found for scalar 
mesons but opposite to that for pseudoscalar mesons. The usual explanation 
for the splitting in the pseudoscalar channel is through the topological charge 
density in the vacuum Isl. 

1 



In lattice studies it is possible to measure separately the non-singlet con- 
tribution which is given by connected correlation C{t) while the flavour singlet 
contribution has an additional disconnected correlation D{t). Previous lattice 
studies have been made of the pseudoscalar mesons 0, pi, 0, 0, 1^ and scalar 
mesons M, nGl. For a discussion including some results for vector and axial 
mesons, see ||ll[ . 

At large t where ground state contributions dominate we have 

C(i) = ce-"i* (1) 

and 

C{t) + D{t) = de-™°* (2) 

where uiq is the flavour singlet mass and mi the flavour non-singlet mass. Now 
if the same meson creation and destruction operators are used for the study of 
both correlations, with quarks degenerate in mass, d and c have the same sign. 
Then by a study oi D/C which is given by 

D/C = (d/c)e("i-™°)* - 1 (3) 

one can explore the mass splitting between flavour singlet and non-singlet. We 
illustrate this behaviour in fig. 1, assuming that only ground state contributions 
occur at all t values. Although it might be thought that d = c, we have shown 
previously jl^l that this is not necessarily the case, and indeed sign changes in 
D/C versus t can be required. So, in summary, the slope (increase/decrease) of 
D/C on a lattice can determine the sign and magnitude of mi — ttiq. 

This has already been explored in detail by us for the cases of pseudoscalar 
mesons M (where mo — mi > ) and scalar mesons |10[ (where niQ — rrii < 0). 
Here we extend the study to vector and axial mesons. 

2 Connected and Disconnected Contributions 

Here we discuss in general what information is available on the connected cor- 
relation C and disconnected correlation D and their relative sign. 

For this section we consider only local operators to create a meson, namely 
•tpATipB, where subscripts refer to the flavour of the quarks. We choose a basis 
of hermitian matrices for F, see Table 1. Here we restrict ourselves to spa- 
tially symmetric operators, a fuller description which takes into account also 
the lattice cubic symmetry is given in |1^ . 

For i > 1 we can use reflection positivity to determine the sign of C. This 
corresponds to using a meson sink -ipB^^iIJA with F^ — 74F^74 = r4F. Then 
the correlation 

c{t) = {i,Ar^B)o{^Br'''iPA)t (4) 
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Figure 1: The ratio of disconnected to connected contributions as given by 
eq-i 



will be positive for i > 1 . We evaluate this correlation in the usual way for the 
connected correlation, combining the Grassmannian fermion fields into prop- 
agators which yields another minus sign, and using the 75 hermitivity of the 
Wilson-Dirac fermion matrix to relate these propagators to a common source. 
We indeed find that C(t) > in our measured correlations for t > 1. 

In the flavour singlet case, when quarks A = B^ there is an additional 
disconnected correlation D(t) to be evaluated. This correlation can be written 
in the form 

I?(i)-iV/r4r5L(0)L*(t) (5) 



where the disconnected loop 



L(t) = TrrAf-i 



(6) 



with M~^ the quark propagator and the sum in the trace is over colour, Dirac 
and spatial indices at time t. 

The factor of rs arises since the Wilson-Dirac fermion matrix M is 75 hermi- 
tian and hence L is real/ imaginary as 75r = r5r75 with r^ = ±1. Now at i = 
we have that L{0)L*{0) > so the disconnected correlation D{0) has sign r4r5. 
If this sign were to be maintained at larger t, then this would give a prediction 
for the sign oi D/C and hence information on the sign of TOi — niQ without any 
lattice evaluation at all. The no-free-lunch is that indeed a sign change in D as 
t increases is possible and it is indeed the goal of this work to explore this on 
the lattice. 

From Table 1, we can deduce that there must be such sign changes in D{t) 
as t increases from 0, since there are two operators available to study both 
pseudoscalar and vector mesons. In each case these two operators have different 
signs of r^r^ and hence one of them must change sign so that they agree on a 
common value of mo — mi at large t where the ground state contribution must 
dominate. This has already been explored for the pseudoscalar case and 
there D(t) for T = 17475 was found to change sign. 
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Table 1: Flavour singlet mesons produced by different operators ipTip and 
the sign factors r^ and r^ as defined in the text whose product determines 
D(Q)/C(t). Here fc is a spatial index and both F and 7 are hermitian. The 
meson quantum numbers are quoted in the continuum limit, for a discussion of 
the appropriate representation of the cubic group on a lattice see Q . We also 
include our observed evidence for sign changes in D{t) from i = which gives 
a prediction of the sign of the non-singlet minus singlet mass shift as shown. 



3 Lattice Methodology 

Here we use dynamical fermion configurations with Nf = 2 from UKQCD |jl3. 
The sea quarks correspond to k = 0.1395 with a tadpole improved clover for- 
malism with Csw — 1.76. 



Local and spatially-fuzzed operators [l3 are used for meson creation (with 
two fuzzed links in a spatially symmetric orientation with 5 iterations of fuzzing 
with coefficient given by 2.5*Straight + Sum of staples). Thus we evaluate a 2 x 2 
matrix of local and fuzzed correlators [O . Mesons created by all independent 
products of gamma matrices are evaluated. Here we restrict our attention to 
the momentum zero sector. 

We measure the disconnected correlations on 252 configurations of size 12'^24 
separated by 20 trajectories. For the evaluation of the disconnected correlators, 
we use stochastic noise sources with variance reduction using the hopping pa- 
rameter expansion [|lO[ . Here we employ more noise sources (96) to explore the 
very small disconnected contributions from some quantum numbers. We use 
sources at every site on the lattice and determine the momentum zero correla- 
tions from them. 

The connected correlator is obtained by explicit inversion from a source (local 
or fuzzed) at the origin for 126 configurations separated by 40 trajectories |13| . 

3.1 Stochastic noise compared to signal 

We measure the zero momentum disconnected loop L{t) on each time-slice for 
each gauge configuration. This ensemble, for each choice of operator F gives us 
the values of the standard deviation CTobs given in Table 2. We also, from our 
96 stochastic samples in each case, have the estimate of the standard deviation 
Cstoch on the mean of these 96 samples coming from the stochastic method. 
We can then deduce the true standard deviation of the gauge time slices from 
fgauge = (o'obs~'^stoch)^''^- This is presented in Table 2. Here the normahsation 
is such that M = I + k . . .. 
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Table 2: Mesons produced by different operators ^/JF?/). The standard deviation 
of the loop operator of eq. o is presented. Here CTstoch is the error estimated from 
the 96 stochastic samples used and this is the used to deconvolute the observed 
spread to give the true standard deviation of the loop (cTgaugc)- 



In an ideal world we would have CTstoch << fgaugo which would imply that 
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Figure 2: The ratio of disconnected to connected contributions for pseudoscalar 
mesons versus time in lattice units. Local and fuzzed operators with F = 75 (+, 
x) and with F = 17475 (diamond, octagon). 



no appreciable error arose from the stochastic methods employed. For the pre- 
viously studied cases, the pseudoscalar with F = 75 and the scalar with F = /, 
we see that the stochastic errors are truly negligible with 96 stochastic samples 
and indeed the 24 samples used before |jiq, m were adequate. For the other 
cases, which have been little studied hitherto, we see that the stochastic errors 
are reasonably small, except for F = 7^ or 74. The latter case has spin exotic 
quantum numbers and is expected to be very poorly determined by our meth- 
ods. For the vector meson, however, one of our goals is to explore the singlet 
mass splitting. Here we see that more stochastic samples (over 1000) would be 
needed with our current method to get the stochastic error significantly smaller 
than the inherent gauge error. We do not have the required computational re- 
sources at present. One small advantage, however, is that we can average over 



the three spin components of the vector which reduces errors by l/v3. Also 
the second operator (F = 274 7fe) which creates a vector meson has a relatively 
smaller stochastic error. We also use fuzzed sources in order to obtain more 
measurements . 
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Figure 3: The ratio of disconnected to connected contributions for vector 
mesons versus time in lattice units. Local and fuzzed operators with F = 7^ (+, 
x) and with F = J747fe (diamond, octagon). 



3.2 Results 

We present in figs. 2 to 4 some of our results for the ratio of the disconnected 
correlator to the connected correlator. The error on the disconnected correlator 
is much larger than that on the connected one. This arises essentially because 
the absolute error on the disconnected correlator stays of the same magnitude 
as t increases, much as is the case for correlations between Wilson loops as used 



in glueball studies. The connected correlator, in contrast, has an approximately 
constant relative error as t increases. For this reason we employ the full data set 
to determine the disconnected correlator: sources at all sites and measurements 
every 20 trajectories. We bin these results to avoid any problem from auto- 
correlations among the gauge configurations separated by only 20 trajectories. 
Even with this approach, the error on the disconnected correlator rises rapidly 
as seen from the figures. A considerably larger number of gauge configurations 
will be needed to explore larger t values. 

Our results for the pseudoscalar case have been presented before and 
are included here for comparison. We do indeed see a consistent slope for the 
two different 7 matrix operators considered with a sign change for T ~ 47475 
as discussed above. The mass difference a{mo — mi) ~ 0.1 is positive here as 
expected since the 77' is heavier than the non-singlet pseudoscalar mesons. 

For the vector mesons, we again have two different 7 matrix combinations 
available. We expect a sign change in one case as discussed above. From fig. 3 
we see that the situation is that the ratio D/C is very small so that detailed 
study is difficult. Moreover the different values observed for fuzzed and local 
operators suggest that excited state contributions are significant, especially for 
the r = 17475 case. The results up to t = 3 do suggest that D/C is negative 
and that a sign change occurs for the F = 7/,. case between t ~ (where it is 
positive) and t — 3. This would imply a(?Tio — rni) > 0. 

For the axial mesons, our results are presented in fig. 4. The /i meson case 
shows a larger signal, consistent with a{mQ — mi) > 0. We find a sign change in 
D{t) from t = Qtot>0 for the hi case, giving a very small signal for D{t)/C{t) 
at small t values which appears to be negative, implying a(mo — mi) > 0. 

We have attempted to make fits to the singlet correlators to extract masses 
in each case, but the rapid increase in errors with increasing t makes these fits 
poorly determined and with results that depend strongly on the minimum and 
maximum t values fitted. A less sensitive way to estimate the allowed magnitude 
of the singlet non-singlet mass splitting is from the values of D{t)/C{t) directly: 
bearing in mind the scenarios illustrated in fig. 1. From the data shown for 
i < 5, we estimate a{mo~m,i) « 0.10(2) for the pseudoscalar case, a{mQ — m,i) « 
0.001(2) for the vector case, a{m,Q — mi) « 0.005(10) for the /i, ai case and 
a{m,Q — m,i) « 0.001(10) for the hi, bi case. Note in particular the hierarchy 
of errors: with smallest errors for the vector case, followed by the axial mesons. 
Note also, as discussed above, that we have information on the sign of a(TOo— mi) 
which is not folded into these error estimates. 

4 Discussion 

Our exploratory study has a scale set by |13] ro/a ~ 3.44 and peudoscalar meson 
to vector meson mass ratio of mpjmy = 0.71. Using the conventional value 
ro = 0.5 fm then gives a^^ = 1.34 GeV while the meson mass ratio implies that 



the sea quarks have masses close to that of a strange quark. Here we are using 
Nf = 2 flavours of degenerate quark in the sea. Since our lattice evaluation is 
for unphysical parameters, we first discuss the experimental spectrum of light 
quark mesons to aid in comparison with our results. 

For a meson considered as made from two valence quarks of either type n (w, 
or d here treated as degenerate) or of type s, we have two observable non-singlet 
states nn and ns. Then the singlet fin state would have a mass degenerate with 
the non-singlet nn if there were no disconnected contributions. Furthermore, 
from the observed non-singlet masses, by assuming an equal splitting in mass 
(or usually in mass squared), one also can deduce the mass which the singlet 
ss state would have if there were no disconnected contributions. Comparison 
of these theoretical mass values with the observed flavour singlet masses, then 
gives information from experiment about the size and nature of disconnected 
contributions. 

Thus there are two observed states in the nonet with isospin zero which can 
have a component of the disconnected contribution and these can mix in terms 
of their quark content. So for flavour singlet sector, we then have contributions 
to the mass squared matrix with quark model content {uu + dd)/\/2 and ss 
(which we label as nn and ss respectively): 

\f2Xns m^^ + Xss J 

Here m corresponds to the mass of the flavour non-singlet eigenstate as 
discussed above and is the contribution to the mass coming from connected 
fermion diagrams while x corresponds to the contribution from disconnected 
fermion diagrams. In the limit of no mixing (all a; = 0, the OZI suppressed 
case), then we have the quenched QCD result that one flavour singlet state is 
degenerate with the isospin one hn meson while the other corresponds to the ss 
meson. 

Using as input r7i„„, nigg and the flavour singlet masses toq and ttiq, the 
three mixing parameters x cannot be fully determined. However if one makes 
some assumption about the mixing parameters x one can deduce the mixing 
pattern - see M for a discussion of this for pseudoscalar mesons. One simple 
case that may be used is to assume that all mixing strengths x are the same for 
mesons of a given J^''^' . Furthermore, if the value of x"^ /{nil^ — nri^^) is small, 
then the off-diagonal contribution to mixing is negligible and the flavour singlet 
mass eigenstates will be of mass squared ?ti^„ -I- 2a; and m^^ -I- x. This gives two 
opportunities to estimate x from the observed spectrum. Also the mass shift of 
the nn flavour singlet (2x) is just the same as in the case of Nf = 2 degenerate 
quarks - which is the case we explore on the lattice here. A caveat applies for 
axial mesons: since charge conjugation is not a good quantum number for the 
sn states, there will be mixing between the J^'~^ = 1++ and 1^ mesons. This 
complicates the mixing scheme further. 



We now present the conclusions of such an analysis based on the experimen- 
tal ||l5| mass values. For the vector mesons, the nn sector gives mo — mi ~ 0.013 
GeV while the ss sector gives 0.016 GeV with a mass squared formalism and 
0.001 GeV with a linear mass formalism. These signs suggest that the splitting 
in the N f = 2 sector would be toq — mi sa 0.01 GeV. 

For the axial mesons, the additional mixing of the ns states only allows an 
analysis if more assumptions are made. Then assuming that the lightest iso- 
singiet state is predominantly nn yields the two flavour result that ttiq — mii « 
0.05, —0.06 GeV for the J^^ = 1^^, 1^ mesons respectively. 

A complication that arises in comparing with experimental meson spectra is 
that of decays. In our lattice studies, since the quark mass is relatively heavy 
(heavier than the strange quark since mp/my = 0.71 while we expect ||7| this 
ratio to be 0.682 for strange quarks) then we have no decay channels open for 
the ground state mesons we study. In contrast, some of these mesons have 
substantial experimental decay widths (150 MeV for the p, over 250 MeV for 
the ai and 140 MeV for the &i). One consequence of this, as has been known 
for a very long time ||l6[, is that the pole in the complex plane corresponding 
to a resonance has an energy whose real part is lower than the quoted value 
which corresponds to a phase shift of 90*^. This mass shift arises from the energy 
dependence of the width and will be more significant for wider resonances. Aside 
from this inherent uncertainty, there may be further dynamical effects arising 
from the back-reaction of the decay channels to the effective propagator. 

Thus we should interpret our results as giving an indication of the strength 
and sign of OZI violating contributions to the light meson spectrum. These 
need not correspond to those observed experimentally because of the above 
issues (namely the more complex mixing schemes allowed for Nf = 3 and the 
decay effects) and also because we would need to extrapolate our lattice results 
to the continuum limit and to more realistic quark masses. In particular, there 
is evidence from the pseudoscalar mesons that the splitting does increase with 
decreasing quarks mass fj, ||, H 0, || . 

We find a hierarchy of singlet non-singlet mass splitting which is large for 
pseudoscalar mesons (0.13(2) GeV), smaller for axial mesons (0.007(13) GeV for 
/i and 0.001(13) GeV for hi) and smallest for vector mesons (0.002(3) GeV). 
This is in agreement with the hierarchy of magnitudes seen experimentally We 
find that in each case the sign of the effect is that the flavour singlet state 
is heavier. This is the sign found experimentally for the pseudoscalar, vector 
and 1"*"+ axial but not for the 1"' axial. The magnitudes we find are smaller 
than the experimental values (except for the hi case where the magnitude is 
comparable), which can come in part from the use of too heavy a light quark 
(as is known to be the case for the pseudoscalar mesons) but also from the 
impact of the large experimental decay widths. 

In this study we have explored from first principles in QCD the OZI rule 
for the meson masses for the case of two degenerate quark flavours. It would 
be interesting to extend this lattice study to OZI rule effects in meson decays. 
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We find that the disconnected contributions (ie OZI violating terms) to the 
masses are indeed small for axial and vector mesons. For the vector and axial 
mesons, we find evidence that the flavour singlet mass is increased compared to 
the non-singlet. This is the opposite of what would be expected in the simplest 
theoretical model: namely mixing with a heavier glueball of the same J^'~^ . An 
understanding of this remains a theoretical challenge. 
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Figure 4: The ratio of disconnected to connected contributions for axial mesons 
(/i and hi types) versus time in lattice units. Local and fuzzed operators for 
the /i meson with F = i"f5"fk (+, x) and for the hi meson with F = i^i'y^lk 
(diamond, octagon). 
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